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ECOLOGY AND THERMAL INACTIVATION OF MICROBES 
IN AND ON INTERPLANETARY SPACE VEHICLE 
COMPONENTS 


INTRODUCTION 

The expression and evaluation of experimental data represent 
one of the more difficult tasks in all areas of science. The 
following manuscript, to be published in Applied Microbiology , 
represents a summary of some of the work done in this laboratory 
relative to statistical techniques of analysis which has relevance 
to studies on the thermal inactivation of bacterial spores. 
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A Comparison of Four Methods Commonly Used 
to Estimate LD50 

J. T. Peeler 

Division of Microbiology 
Food and Drug Administration 
Cincinnati, Ohio 45226 

Biologists are often required to determine the potency of a drug or 
serum, but concentrations are difficult to measure directly by chemical 
methods, and some toxic products may be of unknovm composition. Bio- 
assay is a common technique used to determine the lethal effect of a 
substance and thus compare the results to a known standard. This tech- 
nique requires that test animals (rats, pigs, cats, etc.> be injected 
with several doses of the substance and that the number of survivors at 
each dose be recorded. The number of subjects at each dose is often 
held constant, although this is not necessary. These data are used to 
estimate the mean lethal dose (LD50) , or the dose at which 50 percent 
of the animals survive. 

Ten or more methods are used to estimate the LD50. Finney (5) has 

If 

discussed some of these methods, and he favors the Spearman-Karber and 
moving average angle methods for small sample sizes. He states that 
estimates will, in general, be better by these methods in the face of 
skewed tolerance distributions. 

Two other widely employed methods are the problt analysis (4) and 
the Reed-Muench method (7). The problt method is often used by workers 
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with some statistical training, and the Reed-Muench method is widely 
used by researchers trained in biological sciences. Finney (5) objects 
to ft\e Reed-Muench method because, among other things, the laethod lacks 
an estimate of precision from a single experiment. Brown (2) has given 
an estimate of variance for the Reed-Muench method so that precision can 
now be compared. 

The present study was designed to compare the results of the four 
tests — probit analysis, moving average angle, Reed-Muench, and Spearman- 

tf 

Karber. Three known tolerance distributions were employed (two symmetri- 
cal and one skewed) , and a hundred random samples were generated from 
each. Also, four sample sizes were used, assuming 5, 10, 15, and 20 test 
subjects at a dose. The four tests were used to estimate LD50 and pre- 
cision from the three distributions where the known mean and variance 
could be compared to the estimates. Several sample sizes were used. 
Conclusions were then drawn about methods and sample sizes. 

MATERIALS AND METHODS 

Statistical procedures . The tolerance of Y of any one member of a 
group of experimental animals of a certain species for a specific toxin 
is defined as that dose which would be Just sufficient to cause death; 
with any smaller dose the animal would survive, and with any greater dose 
it would die. The natural logarithm x of y has a probability distribution 
P(x) given by 

P(x) - ' f(t)dt 


— 00 


I 


I 


I 


- 3 - 

which expresses the probability that the "log tolerance" of a randomly 
chosen animal does not exceed x. The probability density function of log 
tolerances, f(t), is generally taken to be the Gaussian or normal density, 
a symmetric function. The mean and median of any symmetric density coin- 
cide; the median of a log tolerance density function is equivalent to the 
logarithm of the median lethal or effective dose, l.e.. In (LD50) . Thus, 

if Y * LD50 and x = In (Y ) , then 
0 0 o 

X 

PCXq) f(t)dt = 0.5 

— 00 

An experimental animal chosen at random and administered the dose Y^ con- 
sequently has a 50% chance of survival. 

The above concept is a reasonable and commonly assumed way of view- 
ing the variation in tolerance to a toxin among the population of all 
animals of a given species. The existence of such a log tolerance dis- 
tribution is assumed in this paper. 

The four methods used in this study were adapted to the computer. 
Calculations for the problt are given by Finney (4) , and computations 
used here were those of Daum et al . (3), where the computer calculates 
the best regression rather than the experimenter choosing values from a 

M 

graph. The Spearman-Karber method was written from calculations given by 
Finney (5) , and the Reed-Muench method was computed as given by Reed-Muench 
(7), with confidence intervals given by Brown (2). Bennett (1) presented 
the modification u.' the moving average angle procedure used here, with 
confidence limits as shown by Harris (6) . 
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Estimates of In (LD50) obtained by Reed-Muench do not lend them- 
selves to any kind of theoretical Investigation as has been made for the 
problt, Spearman-Karber , and moving average methods. Regarding the 
latter two, Finney (5) makes the following recommendations; 

1. When the experimenter knows nothing about the value of the LD50 
of the toxin question, he must necessarily choose a wide dose 
range to be sure of bracketing It. In such cases, the moving 
average method with a span K as large as possible Is preferable 

tl 

to the Spearman-Karber estimate. 

2. If the experimenter Is confident of bracketing the LD50 within a 
small range and d/o Is larger than 1.0, the two methods have 
comparable precision for the same number of animals per dose, n. 

II 

When d/o < 1.0, the precision of the Spearman-Karber method Is 
markedly superior to that of the moving average method. 

II 

3. The experimenter who Intends to employ the Spearman-Karber 
method Is best Instructed to use a few animals (small n^) at 
each of many Oarge m) doses, rather than the same total number 
of animals (N) distributed between fewer doses over the same 
total range. 

Unless n^^ Is so small as to reduce the effective precision of the 
problt estimate on account of a large value of g, the difference between 
the precisions of Spearman-Karber and problt applied to data for which 
each Is suitable will then be negligible. Since nothing Is known about 
the behavior of the precision of the Reed-Muench estimate for varying 
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values of m, and d/o, a comparative study of this estimate with 
those of the other three methods was performed by simulation. The pro- 
cedure followed is described in the following section. The value m is 
the number of doses, n^ the subjects per dose, and d the width between 
doses (i.e. 1 log cycle). 

The simulation procedure . Three tolerance distributions were con- 
sidered ior th -s simulation. One is a mean 4, o *2/^, which is 

skewed. Two normal distributions, which are symmetrical, were chosen 
also, one with mean - - 10.36161, a “ 2.47060 (normal I), and the other 
with mean = “10.36161, o = 6.90780. There are 10 total cases, with sub- 
jects per dose being n = 5, 10, 15, and 20 for normal I and X^(^)* Two 

i 

sample sizes, n^ = 5, 15, are considered for normal II. 

Values for nonsurvivals (r^^) per dose were simulated in the fol.low- 
ing manner. First, a log tolerance probability distribution F and a set 
of m doses were chosen. Then a sample of n^ values of X ■ In y were 
drawn from the chosen F by means of a random number generator in the com- 
puter program. The value r^ Is the number of values X in this sample 
which exceed X^ ■= In Y^ (corresponding to r^ deaths) . This procedure was 
repeated for each dose level. The resulting set of values r^ was then 
used to calculate estimates of In (LD50) by each of the four methods. One 
hundred sets of randomly chosen samples were generated for the 10 cases. 

The normal I case is the Ideal situation for a normally distributed 
In tolerance. The mean of the In tolerance distribution (-10.36161) was 
chosen to be the midpoint of the dosage range, and the o of the distribu- 
tion was chosen so that the range of In dosages was about -10.36161 + 3 o. 



The dosages given in Table I are those values assumed to be examined 
under experimental conditions. In other words, one assumes that the 
experimenter has chosen six dosages and exactly bracketed the true mean 
for the normal I case. 

The normal II case assumes a mean of -10.36161 and o ■ 6.90780, but 
the dosages were chosen so that the true mean falls at the upper **nd of 
the range. This is what might happen when the experimenter tries to 
bracket the In LD50 but is off somewhat. The In tolerance distribution 
is centered at -10.36161 and has a 3 o spread from -31.08501 to 10.36179 

as compared to the examined range of -20.72322 to -9.21032. 

2 

Case three, x (^) * differs from the others in that there exists a 
minimum lethal dosage with an In of zero, and the distribution is skewed 
or asymmetric. The mean and the median values do not coincide. Estimates 
froii! the four methods are thus examined when the true distributions are 
both sjTnmetrical and asymmetrical. The true In (LD50) ■ 4,00000. Dis- 
cussion and results from the simulation procedure are given below. 

RESULTS 

Estimates were obtained for LD50 by the four methods dercribed above. 
The means of 100 random sample estimates are presented in Table 2 for 
each of the 10 cases. Ninety-five-percent confidence limits can be com- 
puted by adding and subtracting the factor given in the table. When 100 
estimates were not computed for a method and case, then the number used 
Is presented. This was often the case for the problt analysis. An inter- 
nal computational rouf'ne was used to test to see if the regression was 
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significant (i.e., to test the null hypothesis that ■ 0 in the 
linear regression). When the regression is not significant, the print- 
out states this result and does not print out estimates. 

The estimates of variance among LD50 values and the coefficients 

of variation expressed as a percent are given in Table 3. These results 

show the expected decrease in variation as the number (n^) of subjects 

per dose Incieases. The coefficient of variation <CV) is satisfactory 

2 

for the nonael I cases, but more than 10% for all x (^) results. This 
means that tcnslderable variation occurs when tlie sampling distribution 
is skewed arj also when the true mean is not centered in the dosage 
range as in normal II. Although the variance decreases with the increas- 
ing subjects per dose, n^ ■ 10 seems Co be a practical number per dose 
to suggest. 

Mean results obtained for the methods can be compared to the known 

It 

means In Table 1. Spearman -Karber tends to underestimate the mean for 
the normal I cases, with the other means being close. The confidence in- 
tervals for Reed-Muench, moving average angle, and problt analyses in 
fact Include the true mean. For the normal II cASaB, Reed-Muench and 
moving average angle methods yield high results, and Spearman-Karber is 
very low In its estimates. Only the problt analysis with larger sample 
size seems adequate. The estimates are closer for the x (^) cases, with 
problt analysis yielding higher estimates than the value of 4.00000 and 
the other three estimates being lower. For n^ ■ 10, only the confidence 
limits of Reed-Muench and problt analysis Include the true mean. The pro- 
bit analysis was found to be the unbiased estimator for all cases. It 
also has the minimum variance for cases where n^^ ^ 10. 
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DISCUSSION 

The results from this simulation study indicate that probit 

analysis is the most desirable method (n^ ^10) . There are internal 

tests of consistency, and the mean LD50 estimates are near the true mean 

(i.e., it is the unbiased minimum variance estimate). In most cases 

probit analysis is also the method with minimum variance. Using estl* 

mates from the other methods for the normal II case might tend to give 

misleading results. The Reed-Muench method does not have any internal 

statistical test, so although it generated an estimate for every sample, 

the LD50 estimates were higher than the true mean more than 90Z of the 

time. All methods performed well for the ideal case where the true mean 

was centered in the dosage range examined. For the normal I , ^ 'v 1 and 

0 

x 2(4) i 'U 0.8. 

Another reason for estimation problem was that for normal II, 

0.3 means that the dosages did not cover the range as they must for the 

II 

Reed-Muench, Spearman-Karber , and moving average angle. For normal I, 
the (-^ 'v 1) dosage range did cover the range from all positive and all 
negative responses, and the plot of estimates was symmetrical about the 
mean, with only a small difference between estimated and true mean. The 
variance estimates (Table 3) indicate that there is little improvement 
above 10 subjects per dose, and this number is suggested as adequate. 

Two restraints were placed on the study from the outset in order to 
control the number of factors to be examined: First, the range of data 

was set as lO'^, lO" etc. and second, the number of doses was re- 

stricted to six. This situation represents the typical type of 
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experiment used by biologists end microbiologists. More dosen would be 
too expensive, as 10 subjects per dose would require 60 for e single 
determination. 

Thus if it is known that the tolerance distribution is symmetrical 
and that the doses contain the true mean, then any one of the four 

II 

methods can be used. Reed-Muench and Spearman-Karber are about the same 
in terms of calculation time, but if the data do not go from all positive 
to all negative, then estimation is poor, as shown by the normal II case. 
Usually, an experimenter will change his choice of doses and run the ex- 
periment again. However, the temptation may be great to use the estimate 
obtained, which may be far from the true mean. Thus, problt analysis 
should be used if possible, and certainly if a computer is available. If 
the experimenter employs fewer than five subjects per dose, however, pro- 
bit analysis should not be used. For cases similar to normal I, the 

II 

Spearman-Karber and Reed-Muench methods give satisfactory results and are 
easy to compute without the aid cf calculators. 
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Table 1. Administered dosages for simulation test! 


Normal I 
dose 

V - -10.36161 
a ' 2.47060 

In (dose) 

Normal II 
dose 

y - 10.36161 
0 - 6.90780 

In (dose) 

dose 

V ■ 4.00000 
0 - 2.82800 
In (dose) 

1 X 10‘^ 

-16.11806 

1 

X 10"^ 

-20.72322 

1 

X 

10 

0.00000 

1 X 10'^ 

-13.81548 

1 

X 10*® 

-18.42064 

1 

X 

10^ 

2.30258 

1 X 10’^ 

-11.51290 

1 

X 10"2 

-16.11806 

1 

X 

2® 

10 

4.60516 

1 X 10~^ 

- 9.21032 

1 

X 10*® 

-13.81548 

1 

X 

10® 

6.90774 

1 X 10‘3 

- 6.907/4 

1 

X 10“- 

-11.51290 

1 

X 

10^ 

9.21032 

1 X 10*2 

- 4.60516 

-1 

X 10*^ 

- 9.21032 

1 

X 

10^ 

11.51290 


a 


True mean lies between these two doses. 



of distribution sampled = -10.36161. 
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Table 3. Variance of LD50's estimated from the 100 random samoles 




